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ABSTRACT 

The globular cluster M15 is unique in its display of star-to-star variations 
in the neutron-capture elements. Comprehensive abundance surveys have been 
previously conducted for handfuls of M15 red giant branch (RGB) and red hor- 
izontal branch (RHB) stars. No attempt has been made to perform a single, 
self-consistent analysis of these stars, which exhibit a wide range in atmospheric 
parameters. In the current effort, a new comparative abundance derivation is pre- 
sented for three RGB and six RHB members of the cluster. The analysis employs 
an updated version of the line transfer code MOOG, which now appropriately 
treats coherent, isotropic scattering. The apparent discrepancy in the previously 
reported values for the metallicity of M15 RGB and RHB stars is addressed and 
a resolute disparity of A{RHB — RGB) ^ 0.1 dex in the iron abundance was 
found. The anti-correlative behavior of the light neutron capture elements (Sr, 
Y, Zr) is clearly demonstrated with both Ba and Eu, standard markers of the 
s- and r-process, respectively. No conclusive detection of Pb was made in the 
RGB targets. Consequently for the M15 cluster, this suggests that the main 
component of the s-process has made a negligible contribution to those elements 
normally dominated by this process in solar system material. Additionally for 
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the M15 sample, a large Eu abundance spread is confirmed, which is comparable 
to that of the halo field at the same metallicity. These abundance results are 
considered in the discussion of the chemical inhomogeneity and nucleosynthetic 
history of M15. 



1. INTRODUCTION 

Detections of multiple main sequences and giant branches in globular clusters (GC; 
e.g. coCen, NGC 2808, and NGC 1851; Bedin et al. 2004(l9D, Piotto et al. 2007(j72D, Han 
et al. 2009 (j^^ ) have challenged the notion that these objects are uniformly mono- metallic 
stellar systems of unique age. In addition to the metallicity variations observed in certain 
clusters, the anomalous abundance behaviors of globulars include star-to-star scatter of light 
element [el/Fe] ratios (for C, N, O, Na, Mg, and Al) in both main sequence and giant 
stars (this is in contrast to the abundance trends of halo field stars; e.g. Carretta et al. 
2009a( 1THj) and Gratton, Sneden, & Carretta 2004 (13U|) )[^ These departures in the relative 
abundances (found in stars of different evolutionary stages) imply that there are multiple 
stellar generations present within the globular cluster and that an initial generation may have 
contributed to the intracluster medium. It is possible that three sources are responsible for 
the aggregate chemical makeup of a globular cluster: a primordial source that generates the 
initial composition of the protocluster cloud, a pollution source that deposits material into the 
ICM from highly-evolved asymptotic branch stars, and a mixing source that is independent 
of the other two and the result of stellar evolution processes. Further discussion of these 
scenarios may be found in, e.g. Bekki et al. (2007) (1T0|) and Carretta et al. (2009b)(|15]). 

On the other hand in the vast majority of globular clusters, minimal scatter in the 
element abundance ratios with Z>20 has been observed. The abundances for the neutron(n- 
) capture elements europium and barium have been measured in several GC's, and only in 
a few, exceptional cases have significantly large intracluster differences in these values been 
seen {e.g. M22; Marino et al. 2009(65)). The predominant mechanism of Eu manufacture is 
the rapid n-capture process (r-process) whereas the primary nucleosynthetic channel for Ba is 
slow n-capture (s-process; additional information pertaining to these production mechanisms 
may be found in e.g. Sneden, Cowan, & Gallino 2008 (1831) ). Consequently, the abundance 
ratio of [Eu/Ba] is used to demonstrate the relative prevalence of the r- or s-process in 



^We adopt the standard spectroscopic notation (Heifer, Wallerstein, & Greenstein 1959(j46)) that for 
elements A and B, [A/B] = logio(NA/NB)* - logio(NA/NB)o. We also employ the definition log e(A) = 
logio(NA/NH) + 12.0. 
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individual stars. In globular clusters with a metallicity of [Fe/H]< —1, a general enhancement 
of [Eu/Ba] ~ +0.4 to +0.6 dex is detected, which indicates that n-capture element production 
has been dominated by the r-process (Gratton, Sneden, & Carretta 2004(|36|) and references 
therein). This in turn is suggestive of explosive nucleosynthetic input from very massive 
stars. 

The very metal-deficient globular cluster M15 (NGC 7078; [Fe/H]~ —2.3) has been 
subject to several abundance investigations including the recent study by Carretta et al. 
(2009a( IT8ll ). They employed both medium-resolution and high-resolution spectra of over 80 
red giant stars to precisely determine the metallicity of this cluster: <[Fe/H]>= —2.314 ± 
0.007. Additionally, they detected variations in the light element abundances (Na and O) 
for stars along the entirety of the Red Giant Branch (RGB). Prior studies have also ob- 
served large scatter in the relative Ba and Eu (intracluster) abundances. With the spectra 
of 17 RGB stars, Sneden et al. (1997(l86|l) found a factor of three spread in both ratios: 
<[Ba/Fe]>= 0.07; a = 0.18 and <[Eu/Fe]>= 0.49; a = 0.200 They were able to exclude 
measurement error as the source for the scatter and determined that the variations were 
correlated: <[Eu/Ba]>= 0.41; cr = 0.11. In a follow-up study of 31 M15 giants by Sneden, 
Pilachowski & Kraft (2000b (i88j)), the scatter in the relative Ba abundance was confirmed: 
<[Ba/Fe]>= 0.12; a = 0.21 (limitations in the spectral coverage did not permit a corre- 
sponding analysis of Eu) . 

The majority of M15 high-resolution abundance analyses have employed yellow-red visi- 
ble spectra to maximize signal-to-noise (stellar flux levels are relatively high for RGB targets 
in this region). In order to precisely derive the neutron capture abundance distribution in 
M15, Sneden et al. (2000a (l85]) ) re-observed three tip giants in the blue visible wavelength 
regime (which contains numerous n-capture spectral transitions). The abundance determi- 
nations of 8 n-capture species (Ba, La, Ce, Nd, Sm, Eu, Gd and Dy) were performed and 
large star-to-star scatter in the all of the [El/Fe] ratios was measured. They also found that 
the three stars exhibited a scaled solar system r-process abundance pattern. Additional 
verification of these abundance results was done by Otsuki et al. (2006(71)) in an analysis 
of six M15 RGB stars (the two studies had one star in common, K462). Consistent with 
Sneden et al. (2000a(l85l)). they detected significant variation in the [Eu/Fe], [La/Fe] and 
[Ba/Fe] ratios. Furthermore, Otsuki et al. found that the ratios of [(Y,Zr)/Eu] show dis- 
tinct anticorrelations with the Eu abundance. Finally employing an alternate stellar sample, 
Preston et al. (2006(73)) examined six red horizontal branch (RHB) stars of MlSj^ For the 



^The anomalously nitrogen-enriched star K969 is omitted; see Appendix A of Sneden et al. 1997(186]) 
^The papers from Sneden et al. (1997), Sneden, Pilachowski, & Kraft (2000b), Sneden et al. (2000a), 
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elements Sr, Y, Zr, Ba and Eu a large (star-to-star) spread in the abundances was measured. 
In essence, all of these investigations have observed considerable chemical inhomogeneity in 
the n-capture elements of the globular cluster M15. 

Two issues are brought to light by the M15 abundance data: the timescale and efficiency 
of mixing in the protocluster environment; and, the nucleosynthetic mechanism(s) responsi- 
ble for n-capture element manufacture. In this globular cluster, large abundance variations 
are seen in the two stellar evolutionary classes as well as in both the light and heavy neutron 
capture species. There is a definitive enhancement of r-process elements found in some stars 
of M15 [e.g. K462), yet not exhibited in others [e.g. B584). Taking into consideration the 
entirety of the M15 n-capture results, these data hint at the existence of a nucleosynthetic 
mechanism different from the classical r- and s-processes. Evidence of such a scenario (with 
multiple production pathways) may be also found in halo field stars of similar metallicity 
such as CS 22892-052 (Sneden et al. 2003(j8lD) and HD 122563 (Honda et al. 2006(I17D), 
which have displayed similar abundance variations. Indeed, several models have advanced 
the notion of more than one r-process formation scenario {e.g. Wasserburg & Qian 2000, 
2002 ([SSI EH), Thielemann et al. 2001(1101), and Kratz et al. 2007(1S51)). 

To further understand the implications of the M15 results, the spectra from the three 
RGB stars of Sneden et al. (2000a) dHS} and the six RHB stars of Preston et al. (2006) (73) 
are re-analyzed. A single consistent methodology for the analysis is employed and an ex- 
pansive set of recently-determined oscillator strengths is utihzed {e.g. Lawler et al. 2009 (1^ . 
Sneden et al. 2009 (1871) . and references therein). As the pre- and post-He core fiash giants 
are examined, the relative invariance of abundance distributions will be ascertained for r— 
and s-process species. In consideration of the M15 investigations cited above, a few data 
anomalies have come to light. The two main issues to be resolved include: large discrepan- 
cies in the log e{El) values between the studies of Sneden et al. (2000a) (1551) and Otsuki et 
al. (2006) (I7T1) : and, the significant disparity in the derived metallicity for the M15 cluster 
between Preston et al. 2006 {<[Fe/'H]>RHB= —2.63) and the canonically accepted value of 
<[Fe/H]>= -2.3 {e.g. Carretta et al. 2009a(ll8l), Sneden et al. 2000a(l85D). It is suggested 
that these differences are mostly due to selection of atomic data, model atmosphere, and 
treatment of scattering. 



and Preston et al. (2006) are from collaborators affiliated with institutions in both California and Texas. 
Hereafter, these papers and other associated publications will be referred to as CTG. 



- 5 - 



2. OBSERVATIONAL DATA 

For the three RGB stars of the M15 cluster, the re-analysis of two sets of high resolution 
spectra was performed: the first from Sneden et al. (1997 (1861) ) with approximate wavelength 
coverage region of 5400A < A < 68OOA and the second from Sneden et al. (2000a (1551) ) with a 
wavelength domain of 3600A < A < 5200A. All spectral observations were acquired with the 
High Resolution Echelle Spectrometer (HIRES; Vogt et al. 1994({9l)) at the Keck I 10.0-m 
telescope (with a spectral resolving power of R = A/ A A ^ 45000). The signal-to-noise (S/N) 
range of the data varied from 30 < S/N < 150 for the shorter wavelength spectra to 100 
< S/N < 150 for the longer wavelength spectra (the S/N value generally increased with 
wavelength). The three giants, K341, K462, and K58iQ were selected from the larger stellar 
sample of Sneden et al. (1997) fl86p due to relative brightness, rough equivalence of model 
atmospheric parameters, and extreme spread in associated Ba and Eu abundances. 

Re-examination of the high resolution spectra of six RHB stars from the study of Preston 
et al. (2006) C73D was also done. The observations were taken at the Magellan Clay 6.5-m 
telescope of the Las Campanas Observatory with the Magellan Inamori Kyocera Echelle 
(MIKE) spectrograph (Bernstein et al. 2003) (pi). The data had a resolution of R~ 40000 
and the S/N values ranged from S/N ~25 at 3600 A to S/N ~120 at 7200 A (note that 
almost complete spectral coverage was obtained in the region 3600 A < A < 72OOA). The SIX 
RHB targets were chosen from the photometric catalog of Buonanno et al. (1983) (ll3]) and 
accordingly signified as B009, B028, B224, B262, B412 and B584. It should be pointed out 
that these stars have significantly lower temperatures than other HB members (and thus, 
match up favorably with the RGB). 

Figure [l] features the color-magnitude diagram (CMD) for the M15 globular cluster 
with a plot of the V versus {V — K) magnitudes. The V magnitudes for the RGB stars 
are taken from the prelimiary results of Cudworth (2011) and verified against the data 
from Cudworth (1976 (l25p ). Alternatively, the RHB V magnitude values are obtained from 
Buonanno et al. (1983 (IT3j) ). The K magnitudes for all M15 targets are taken from the Two 
Micron All Sky Survey (2MASS; Skrutskie et al. 2006(|HI])). Cluster members with both B 
and V measurements from Buonanno et al. are displayed in the plot (denoted by the black 
dots) and the stars of the current study are indicated by large, red circles. Note that the 
identifications of RGB and RHB members are based upon stellar atmospheric parameters as 
well as the findings from Sneden et al. (1997(186|)) and Preston et al. (2006(73); please consult 
those references for additional details). Also in Figure [l| two isochrone determinations are 
overlayed upon the photometric data: Marigo et al. (2008(64); with the age parameter set to 



■*The Kustner (1921) (58) designations are employed throughout the text. 
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12.5 Gyrs and a metallicity of [M/H]=-2.2; shown in green) and Dotter et al. (2008 (j27|) : with 
the age parameter set to 12.5 Gyrs and a metalhcity of [M/H]=-2.5; shown in blue). These 
are the best-fit isochrones to the general characteristics ascribed to M15 and no preference 
is given to either source. 

Additional observational details of the aforementioned data samples may be found in 
the original Sneden et al. (1997,2000a) (jHSJ ES} and Preston et al. (2006) (jT^j) pubhcations. 
These papers also contain descriptions of the data reduction procedures, in which standard 
IRAF0 tasks were used for extraction of multi-order wavelength-calibrated spectra from the 
raw data frames, and specialized software (SPECTRE; Fitzpatrick & Sneden 1987 (130]) ) was 
employed for continuum normalization and cosmic ray elimination. 

Figure [2] features a comparison of the spectra of all M15 targets. Displayed in this 
plot is a small wavelength interval 4121 — 4133 A, which highlights the important n-capture 
transitions La II at 4123.22 A and Eu II at 4129.72 A. The spectra are arranged in decreasing 
Teff from the top to the bottom of the figure. As shown, the combined effects of T^s and 
log g influence the apparent line strength, and accordingly, transitions which are saturated 
in the RGB spectra completely disappear in the warmer RHB spectra. 

3. METHODOLOGY AND MODEL DETERMINATION 

Several measures were implemented in order to improve and extend the efforts of Sneden 
et al. (1997, 2000a) (IHSl ESI and Preston et al. (2006(j7SI)). Fii'st, the modification of the line 
analysis program MOOG was performed to accurately ascertain the relative contributions to 
the continuum opacity (especially necessary for the bluer wavelength regions and the cool, 
metal-poor RGB targets). Second, the employment of an alternative grid of models was done 
to obtain an internally consistent set of stellar atmospheric parameters for the total M15 
sample. Third, the utilization of the most up-to-date experimentally and semi-empirically- 
derived transition probability data was done to determine the abundances from multiple 
species. 

3.1. Atomic Data 

Special effort was made to employ the most recent laboratory measurements of oscillator 
strengths. When applicable, the inclusion of hyperfine and isotopic structure was done for the 



^IRAF is distributed by NOAO, which is operated by AURA, under cooperative agreement with the NSF. 
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derivation of abundances. Tables 2 and 3 list the various literature sources for the transition 
probability data. Some species deserve special comment. The Fe transition probability values 
are taken from the critical compilation of Fuhr & Wiese (2006 (l33|l ; note for neutral Fe, the 
authors heavily weigh the laboratory data from O'Brian et al. 1991 (l69|l ). No up-to-date 
laboratory work has been done for Sr, and so, the adopted gf-values are from the semi- 
empirical study by Brage et al. (1998(12); these values are in good agreement with those 
derived empirically by Gratton & Sneden 1994(j37])). Similarly, the most recent laboratory 
effort for Y was by Hannaford et al. (1982(H3il). Yet these transition probabilities appear to 
be robust, yielding small line-to-line scatter. 

A particular emphasis of the current work is the n-capture element abundances, for 
which a wealth of new transition probability data have become recently available. Corre- 
spondingly, the extensive sets of rare earth 5f/-values from the Wisconsin Atomic Physics 
Group were adopted (Sneden et al. 2009 (l87|l . Lawler et al. 2009 (162|) . and references therein). 
These data when applied to the solar spectrum yield photospheric abundances that are in 
excellent agreement with meteoritic abundances. For neutron capture elements not stud- 
ied by the Wisconsin group (which include Ba, Pr, Yb, Os, Ir, and Th), alternate literature 
references were employed (and these are accordingly given in the two aforementioned tables). 



3.2. Consideration of Isotropic, Coherent Scattering 

In the original version of the line transfer code MOOG (Sneden 1973 (1821) ). local ther- 
modynamic equilibrium (LTE) was assumed and hence, scattering was treated as pure ab- 
sorption. Accordingly, the source function, S\, was set equal to the Planck function, Bx(T), 
which is an adequate assumption for metal-rich stars in all wavelength regions. However 
for the extremely metal-deficient, cool M15 giants, the dominant source of opacity switches 
from B.~BF to Rayleigh scattering in the blue visible and ultraviolet wavelength domain 
(A < 4500A). It was then necessary to modify the MOOG program as the LTE approxima- 
tion was no longer sufficient (this has also been remarked upon by other abundance surveys, 
e.g. Johnson 2002(150]) and Cayrel et al. 2004(120])). 

The classical assumptions of one-dimensionality and plane-parallel geometry continue 
to be employed in the code. Now with the inclusion of isotropic, coherent scattering, the 
framework for solution of the radiative transfer equation (RTE) shifts from an initial value to 
a boundary value problem. The source function then assumes the fornJjof S = {1 — e)J + eB 



^To re-state, the equation terms are defined as follows: S is the source function, e is the thermal coupling 
parameter, J is the mean intensity, and B is the Planck function. 
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and the description of line transfer becomes an integro-differential equation. The chosen 
methodology for the solution of the RTE (and the determination of mean intensity) is the 
approach of short characteristics that incorporates aspects of an accelerated convergence 
scheme. In essence, the short characteristics technique employs a tensor product grid in 
which the interpolation of intensity values occurs at selected grid points. The prescription 
generally followed was that from Koesterke et al. (2002 (152(1 and references therein). The 
Appendix provides more detail with regard to the MOOG program alterations. 

Prior to these modifications, for a low temperature and low metallicity star {e.g. a RGB 
target), the ultraviolet and blue visible spectral transitions reported aberrantly high abun- 
dances in comparison to those abundances found from redder lines. With the implementation 
of the revised code, better line-to-line agreement is found and accordingly, the majority of 
the abundance trend with wavelength is eliminated for these types of stars. Note for the 
RHB targets, minimal changes are seen in abundances with the employment of the modified 
MOOG program (as the dominant source of opacity for these relatively warm stars is always 
H over the spectral region of interest). 

3.3. Atmospheric Parameter Determination 

To obtain preliminary estimates of Teg and log g for the M15 stars, photometric data 
from the aforementioned sources (Cudworth 2011; Buonanno et al. 1983(13); 2MASS) were 
employed as well as those data from Yanny et al. 1994(40). To transform the color, the color- 
Tcfj relations of Alonso et al. (1999) (jT]) were used in conjunction with the distance modulus 
((m — M)o = 15.25) and reddening {E{B — V^) = 0.10) determinations from Kraft & Ivans 
(2003 (jMI))- Note that an additional intrinsic uncertainty of about 0.1 dex in log g remains 
among luminous RGB stars owing to stochastic mass loss of order 0.1 dex. Consequently, 
initial masses of 0.8 Mq and O.GMq were assumed for RGB and RHB stars respectively. The 
photometric V and (V-K) values as well as the photometrically- and spectroscopically-derived 
stellar atmospheric parameters are collected in Table [T} 

With the use of the spectroscopic data analysis program SPECTRE (Fitzpatrick & 
Sneden 1987(130])), the equivalents widths (EW) of transitions from the elements Ti I/II, Cr 
I/II, and Fe I/II were measured in the wavelength range 3800-6850 A. The preliminary Tg^ 
values were adjusted to achieve zero slope in plots of Fe abundance (log {epei)) as a function 
of excitation potential (x) and wavelength (A). The initial values of log g were tuned to 
minimize the disagreement between the neutral and ionized species abundances of Ti, Cr, 
and Fe (particular attention was paid to the Fe data). Lastly, the microturbulent velocities 
Vt were set as to reduce any dependence on abundance as a function of EW. Final values of 
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TefT, log g, ft, and metallicity [Fe/H] are listed in Table 1, as well as those values previously 
derived by Sneden et al. (2000a(85)) for the RGB stars and Preston et al. (2006) ([73]) for the 
RHB stars. 



3.4. Selection of Model Type 

To conduct a standard abundance analysis under the fundamental assumptions of one- 
dimensionality and local thermodynamic equilibrium (LTE), two grids of model atmospheres 
are generally employed: Kurucz-Castelli (Castelli &; Kurucz 2003 (il9i): Kurucz 2005( 1^7]) ) and 
MARCS (Gustafsson et al. 2OO8(jlI]))0 The model selection criteria were as follows: the 
reconciliation of the metallicity discrepancy between the RGB and RHB stars of M15, the 
derivation of (spectroscopic-based) atmospheric parameters in reasonable agreement with 
those found via photometry, and the attainment of ionization balance between the Fe I and 
Fe II transitions. For the RHB targets, interpolated models from the Kurucz-Castelli and 
MARCS grids were comparable and yielded extremely similar abundance results. However, 
there are a few notable differences between the two model types for the RGB stars with 
regard to the Pgas and Peiectron content pressures. Though beyond the scope of the current 
effort, it would be of considerable interest to examine in detail the exact departures between 
the Kurucz-Castelli and MARCS grids. To best achieve the aforementioned goals for the 
M15 data set, MARCS models were accordingly chosen. 



3.5. Persistent Metallicity Disagreement between RGB and RHB stars 

For the RHB stars, the presently-derived metallicties differ slightly from those of Preston 
et al. (2006): <[Fe7/H]> = -2.69 (a change of A = -0.03) and <[Fe///H]> = -2.64 (a change 
of A = —0.04). However for the RGB stars, the [Fe/H] results of the current study do vary 
significantly from those of Sneden et al. (2000): <[Fe//H]> = -2.56 (a downwards revision of 
A = —0.26) and <[Fe/7/H]> = -2.53 (a downwards revision of A = —0.28). The remaining 
metallicity discrepancy between the RGB and RHB stars is as follows: A{RGB — RHB)pei 
= 0.13 and A{RGB - RHB)Feii = 0.11. Even with the employment of MARCS models and 
the incorporation of Rayleigh scattering (not done in previous efforts), the offset persists. 
Repeated exercises with variations in the T^s, log g, and Vt values showed that this metal- 
licity disagreement in all likelihood cannot be attributed to differences in these atmospheric 



^Kurucz models are available through the website: 



http://kurucz.harvard.edu/ 



be downloaded via the website: 'http : / / marcs . astro .uu. se / 



and MARCS models can 
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parameters. As a further check, the derivation of [Fe/^/7/H] values was performed with a 
hst of transitions satisfactorily measurable in both RGB and RHB spectra. No reduction in 
the metallicty disagreement was seen as the offsets were found to be: A{RGB — RHB)pei 
= 0.11 (with 45 candidate Fe/ lines) and A{RGB - RHB)Feii = 0.14 (with 3 candidate Fe 
II// lines). 

The data from the M15 RGB and RHB stars originate from different telescope/instrument 
set-ups. Additionally, somewhat different data reduction procedures were employed for the 
two samples. Possible contributors to the iron abundance offset could be the lack of con- 
sideration of spherical symmetry in the line transfer computations and the generation of 
sufficiently representative stellar atmospheric models for these highly evolved stars (which 
exist at the very tip of the giant branch and/or have undergone He-core flash). Indeed, it 
is difficult to posit a single, clear-cut explanation for the disparity in the RGB and RHB 
[Fe/H] values. In an analysis of the globular cluster M92, King et al. (1998(51j)) derived an 
average abundance ratio of <[Fe/H]>= -2.52 for six subgiant stars, a factor of two lower 
than the <[Fe/H]> value measured in the red giant stars. Similarly, Korn et al. (2007) (53) 
surveyed turn off (TO) and RGB stars of NGC 6397 and found a metallicity offset of about 
0.15 dex (with the TO stars reporting consistently lower values of [Fe/H]). They argued that 
the TO stars were afflicted by gravitational settling and other mixing processes and as a 
result, the Fe abundances of giant stars were likely to be nearer to the true value. While the 
TO stars do have T^s values close to that of the M15 RHB stars, they have surface gravities 
and lifetimes that are considerably larger. Accordingly, it is not clear if the offset in M15 
has a physical explanation similar to that proposed in the case of NGC 6397. 



4. ABUNDANCE RESULTS 

For the extraction of abundances, the two filters of line strength and contaminant pres- 
ence were used to assemble an effective line list. Abundance derivations for the majority 
of elements employed the technique of synthetic spectrum line profile fitting (accomplished 
with the updated MOOG code as described in §3.2). For a small group of elements (those 
whose associated spectral features lack both hyperfine and isotopic structure), the simplified 
approach of EW measurement was used (completed with both the MOOG code and the 
SPECTRE program; Fitzpatrick & Sneden 1987 (1BU|) ). Presented in Tables |2] and |3] are the 
log e abundance values for the individual transitions detected in the M15 RGB and RHB 
stars respectively. These tables also list the relevant line parameters as well as the associated 
literature references for the (7 /-values employed. 

In addition to the line-to-line scatter, errors in the abundance results may arise due 



- 11 - 



to uncertainties in the model atmospheric parameters. To quantify these errors in the M15 
data set, a RGB target, K462, is first selected. If alterations of ATg// = ±100 K are applied, 
then the abundances of neutral species change by approximately A[X7/H] ~ ±0.15 whereas 
the abundances of singly-ionized species change by about A[X///H] ~ ±0.04. Variations 
of Alog g = ±0.20 yield A[X//H] ~ ±0.04 in neutral species abundances and A[X///H] 
~ ±0.05 in singly-ionized species abundances. Changes in the microturbulent velocity on 
the order of Aft = ±0.20 result in abundance variations of A[X//H] ~ ±0.08 and A[X///H] 
~ ±0.04. The exact same procedure is then repeated for the RHB star, B009. Modifications 
of the temperature by ATe// = ±100 K lead to abundance changes of A[X//H] ^ ±0.09 in 
neutral species and A[X///H] ±0.02 in singly-ionized species. Alterations of the surface 
gravity by Alog g = ±0.20 engender variations of A[X//H] ^ ±0.01 and A[X///H] ^ ±0.07. 
Finally, variations of Aft = ±0.20 produce abundance changes of A[X//H] ^ ±0.08 and 
A[X7//H] ~ ±0.04. 

To discuss the abundance results in the following subsections, the elements are divided 
into four groupings: light (Z = 8; 11 < Z < 21; Figurejs]), iron-peak (22 < Z < 28; Figurejl]), 
light /intermediate n-capture (29 < Z < 59; Figure [s]), and heavy/other (60 < Z < 70; 
Z = 72,76,77,90; Figure [6]). The measurement of abundances was completed for a total 
of 40 species. Note that for the elements Sc, Ti, V, and Cr, abundance determinations 
were possible for both the neutral and first-ionized species. In light of Saha-Boltzmannn 
calculations for these elements, greater weight is given to the singly-ionized abundances (i.e. 
for the stars of the M15 data set, only a small fraction of these elements predominantly 
reside in the neutral state). 

Figures |3| |4| [5} and |6] exhibit the abundance ratios for the M15 sample in the form 
of quartile box plots. These plots show the interquartile range, the median, and the mini- 
mum/maximum of the data. Outliers, points which have a value greater than 1.25 times the 
median, are also indicated. For all of the figures, RGB abundances are signified in red while 
RHB abundances are denoted in blue. Note that the plots depict the abundance results in 
the [Elem/H] form in order to preclude erroneous comparisons of the RGB and RHB data, 
which would arise from the iron abundance offset between the two groups. 

Table |4] contains the <[Elem/Fe]> values for elements analyzed in the M15 sample along 
with the line-to-line scatter (given in the form of standard deviations), and the number 
of lines employed. The subsequent discussion will generally refer to these table data and 
as is customary, present the relative element abundances with associated a values. The 
reference solar photospheric abundances {without non-LTE correction) are largely taken from 
the three-dimensional analyses of Asplund et al. (2005, 2009(0 [6])) and Grevesse et al. 
(2010(l39|)). However, the photospheric values for some of the n-capture elements are obtained 
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from other investigations {^e.g. Lawler et al. 2009 (162]) ). Table [s] lists all of the chosen 
log numbers. Note that in the derivation of the relative element abundance ratios [X/Fe], 
<[Fe//H]> are employed for the neutral species transitions while <[Fe///H]> are used for the 
singly-ionized lines. This is done in order to minimize ionization equilibrium uncertainties 
as described in detail by Kraft & Ivans (2003) CS^I). 

4.1. General Abundance Trends 

Within the RGB sample, the neutron capture element abundances of K462 are consis- 
tently the largest whereas those of K583 are the smallest. The two RHB stars, B009 and 
B262, exhibit abundance trends similar to those of the RGB objects. The expected anti- 
correlations in the proton capture elements {^e.g. Na-0 and Mg-Al) are seen. The greatest 
abundance variation with regard to the entire M15 data set is found for the neutron cap- 
ture elements. Indeed, the star-to-star spread for the majority of n-capture abundances is 
demonstrable for all M15 targets and is not likely due to internal errors. 

Inspection of Table |4] data indicates that RHB stars generally have higher r-process 
element abundances than RGB stars (on average A[Elem/Fe]/jj7B-_RGB ~ 0.3 dex). A size- 
able portion of the discrepancy is attributable to the difference in the iron abundances as 
<[Fe/H]>/j//B is approximately 0.12 dex lower than <[Fe/H]>RGB- The remaining offset is 
most likely a consequence of the small number of targets coupled with a serious selection 
effect. The original sample of RHB stars from Preston et al. (2006 (1731) ) was chosen as a 
random set of objects with colors and magnitudes representative of the red end of the HB. 
These objects were selected without prior knowledge of the heavy element abundances. On 
the other hand, the three RGB stars from Sneden et al. (2000a (l85i) ) were particularly chosen 
as representing the highest and lowest abundances of the r-process as predetermined in the 
17 star sample of Sneden et al. (1997 (l86j) ). 

4.2. Light Element Abundances 

The finalized set of light element abundances include: C, O, Mg, Al, Si, Ca, and Sc////. 
In general, an enhancement of these element abundances relative to solar is seen in the entire 
M15 data set. 

An underabundance of carbon was found in one RHB and three RGB targets of M15 
based on the measurement of CH spectral features. As the forbidden O I lines were detectable 
only in RGB stars, the average abundance ratio for M15 is <[0/Fe]>RGB= +0.75. This 
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value is substantially larger than that found by Sneden et al. (1997) (l86|) . A portion of the 
discrepancy is due to the approximate 0.15 dex difference in the <[Fe7-/H]> values between 
the two investigations. The remainder of departure may be attributed to the adoption of 
different solar photospheric oxygen values: the current study employs log e{0)Q = 8.71 (Scott 
et al. 2009(1791) while Sneden et al. uses log e{0)Q = 8.93 (Anders & Grevesse 1989^). 

For the determination of the sodium abundance, the current study relies solely upon the 
Di resonance transitions. Table |4] lists the spuriously large spreads in the Na abundance for 
both the RGB and RHB groups. The Na Di lines are affected by the non-LTE phenomenon 
of resonance scattering (Asplund 2005b (j4|; Andrievsky et al. 2007([3])), which MOOG does 
not take into account. Also, Sneden et al. (2000b (l88|) ) made note of the relative strength 
and line profile distortions associated with these transitions and chose to discard the [Na/Fe] 
values for stars with Tes > 5000 K. Consequently, the sodium results from the current study 
are given little weight and are not plotted in Figure [3j 

Aluminum is remarkable in its discordance: <[Al//Fe/]>= 0.37 for one RHB and three 
RGB targets whereas <[Al//Fe/]>= —0.43 for five RHB stars. Now, the relative aluminum 
abundances for the RHB stars match well with the values found by Preston et al. (2006( 17311 ). 
Similarly, the Al abundances from the current analysis agree favorably with the RGB data 
from Sneden et al. (1997 (l86j) ). Though relatively strong transitions are employed in the 
abundance derivation, the convergence upon two distinct [Al/Fe] values is nontrivial and 
could merit further exploration. 

A decidedly consistent Ca abundance ratio is found for the RGB sample: < [Cai/Fci] >rgb= 
0.29; and, also for the RHB sample: <[Ca//Fe/]>ij//B= 0.53. After consideration of the iron 
abundance offset, the RHB stars still report slightly higher calcium abundances than the 
RGB stars. Overall, a distinct overabundance of Ca relative to solar is present in the M15 
cluster. Note that the Sc/ abundance determination was done for only one M15 star (and 
gives a rather aberrant result compared to the Sc// abundance data from the other M15 
targets) . 

4.3. Iron Peak Element Abundances 

The list of finalized Fe-peak element abundances consists of Ti////, Vi/n, Ctj/u, Mn, Co, 
and Ni. Due to RGB spectral crowding issues, derivations of [Ni//Fe/] ratios are performed 
only for RHB stars. 

Achievement of ionization equilibrium did not occur for any of the perspective species: 
Ti////, V////, or Cr ////. In consideration of the entire M15 data set, the best agreement 
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between neutral and singly-ionized species arises for titanium, with all <[Ti/Fe]> ratios being 
supersolar. The V// relative abundances compare well with one another for the RGB and 
RHB targets (comparison for V/ is not possible as there are no RHB data for this species). 
Both of the RGB and RHB <[Cr//Fe/]> ratios are underabundant with respect to solar and 
the neutral chromium values match almost exactly with one another (after accounting for 
the [Fe/H] offset). On the other hand, the worst agreement is found for Ctj/jj in RHB stars 
with A(// - /) = 0.47. 

Subsolar values with minimal scatter were found for the <[Mn7'/Fe7']> ratios in both the 
RGB and RHB stellar groups. However in comparison to RGB stars, manganese appears to 
be substantially more deficient in RHB targets. The discrepancy may be attributed to both 
the RGB/RHB iron abundance disparity as well as the employment of the Mn^ resonance 
transition at 4034.5 A for the RHB abundance determination. In particular, Sobeck et al. 
(2011 (1891) ) have demonstrated that the manganese resonance triplet (4030.7, 4033.1, and 
4034.5 A) fails to be a reliable indicator of abundance. Consequently, the RHB abundance 
results for Mn/ are given little weight and are not plotted in Figure |4} 

4.4. Light and Intermediate n-Capture Element Abundances 

Finalized abundances for the light and intermediate n-capture elements include Cu, 
Zn, Sr, Y, Zr, Ba, La, Ce, and Pr. In general, the RGB element abundance ratios are 
slightly deficient with respect to the RHB values. Also, enhancement with respect to solar 
is consistently seen in all M15 targets for the elements Ce and Pr. 

An extremely underabundant copper abundance relative to solar was found in the RGB 
stars: <[Cu//Fe/]>ii;GB= —0.91. A similar derivation could not take place in the RGB 
targets as the Cu/ transitions were too weak. A large divergence between RGB and RHB 
stellar abundances exists for zinc. Detection of the Zn transitions was possible in only one 
RHB target, which could perhaps account for some of the discrepancy. 

For the entire M15 data set, Y// exhibits lower relative abundance ratios in juxtaposition 
to both Sin and Zr//. With regard to the three average elemental abundances (of Sr, Y, Zr), 
moderate departures between the RGB and RHB groups are seen. Also, a large variation in 
the <[Sr///Fe//]> ratio was found for the members of the RGB group. 

Though different sets of lines are employed, the RGB and RHB <[Ba///Fe//]> ratios 
are consistent with one another. A portion of the RHB abundance variation is due to the 
exclusive use of the resonance transitions in the determination (two lowest temperature 
RHB stars report quite high a values; these strong lines could not be exploited in the RGB 
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analysis). Notably for this element group, the greatest star-to-star abundance scatter was 
found for lanthanum: Argb = 0.46 and Arhb = 0.61 (excluding the one RHB outlier). The 
relative cerium abundances also exhibit a wide spread in the RGB sample. 

4.5. Heavy n- Capture Element Abundances 

The list of finalized <[E1/Fe]> ratios for the heavy n-capture elements is as follows: 
Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Hf, Os, Ir, Pb and Th. All of these element 
abundance ratios are enriched with regard to the solar values. As shown in Figure |6} a larger 
abundance spread is found for this group in comparison to the other element groups. 

Note that as Tefj increases, the strength of the heavy element transitions rapidly de- 
creases and as a consequence, the use of these lines for abundance determinations in the 
warmest stars becomes unfeasible. It was possible to obtain robust abundances for Nd, 
Sm, Tb, and Tm in a single RHB target. On the other hand, abundance extractions for 
the species Os and Ir were done only in RHB stars (measurements of these element transi- 
tions were attainable as less spectral crowding occurs in these stars). Nonetheless, minimal 
line-to-line scatter is seen for the bulk of RGB and RHB n-capture abundances. 

A rigorous determination of the europium relative abundance was performed for all M15 
stars: <[Eu///Fe//]>RGB= 0.53 and <\Ej\xn /¥eii\> rhb= 0.88. Despite the iron abundance 
offset, the largest departure between the two stellar groups is found for the element Ho//. 
Further, the greatest star-to-star scatter in the heavy n-capture elements is seen for the 
[Yb///Fe//] ratio: /\rgb = 0.55 and Arhb = 0.58. 

4.6. Comparison with Previous CTG efforts and Otsuki et al. (2006) 

These new abundance results are now compared to those from the four prior CTG 
publications. For the majority of elements, the current data are in accord with the findings 
of Sneden et al. (1997, 2000a(lMl ES}) and Preston et al. (2006(73)). In this effort, abundance 
derivations are performed for 13 new species: Sc/, V//, Cu/, Pin, Tb//, Ho//, Er//, Tm//, 
Yb//, Hf//, Os/, It I, and Pb//. For elements re-analyzed in the current study, the abundance 
data have been improved with the use of higher quality atomic data, additional transitions, 
and a revised version of the MOOG program. A few large discrepancies in the [El/Fe] ratios 
do occur between the current study and the previous M15 efforts. These departures can be 
attributed to the employment of different [Fe/H] and solar photospheric values as well as the 
updated MOOG code. Accordingly, the results from the current analaysis supersede those 
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from the earlier CTG papers. 

As in Sneden et al. 1997(jMD, ^^e abundance behavior of the proton capture elements 
appears to be decoupled from that of the neutron capture elements. Notably for M15, 
significant spread in the abundances was confirmed for both Ba and Eu. The scatter of 
Aloge(-Ba) = 0.48 and Aloge(-E''u) = 0.90 from the current effort is in line with that of 
Aloge(5a) = 0.60 and Aloge{Eu) = 0.73 from Sneden et al. (1997)(!M[). 

Comparison of the findings from the current study to those from Otsuki et al. (2006) (71) 
has also been done and will be limited to the only star that the two investigations have 
in common, K462. Due to differences in the <[Fe/H]> values, the log (e) data of the 
two analyses are compared. The model atmospheric parameters for K462 differ somewhat 
between the current effort (Teg/log g/vt = 4400/0.30/2.00) and Otsuki et al. (Tes/log g/vt 
= 4225/0.50/2.25). However, the agreement in the abundances for the elements Y, Zr, 
Ba, La and Eu is rather good between the two studies, with the exact differences ranging: 
0.01 <| A{Otsuki— Current) \< 0.16. The largest disparity occurs for Sr, with both analyses 
employing the resonance transitions. As mentioned previously, these lines are not the most 
rigorous probes of abundance. 



4.7. General Relationship of Ba, La, and Eu Abundances 

Sneden et al. (1997)( 186|) claimed to have found a binary distribution in a plot of [Ba/Fe] 
versus [Eu/Fe], with 8 stars exhibiting relative Ba and Eu abundances approximately 0.35 
dex smaller than the remainder of the M15 data set. To re-examine their assertion. Figure [7] 
is generated, which plots [(Ba, La)/H] as a function of [Eu/H] for the entire data sample 
of the current study. It also displays the re-derived/re-scaled Ba, La, and Eu abundances 
for all of the giants from the Sneden et al. (1997) (l86|) publication. No decisive offset is 
evident in either panel of Figure [7} For completeness, the equivalent width data from Otsuki 
et al. (2006) (jTTj) were also re-analyzed and the abundances were re-determined. Again, no 
bifurcation was detected in the Ba and Eu dataEl 



®To avoid duplication, the stars from Ostuki et al. are not plotted as they are a subset of the original 
sample from the Sneden et al. (1997 (186] )) study. 
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5. DATA INTERPRETATION AND ANALYSIS 

A significant amount of r-process enricliment lias occurred in tlie M15 globular cluster. 
Figure [sjplots the average log e values of the n-capture elements (with 39 < Z < 70) for three 
RGB stars (K341, K462, K583; signified by red symbols) and three RHB stars (B009, B224, 
B262; denoted by blue symbols)]^ The solid black line in this figure indicates the scaled, 
solar r-process prediction as computed by Sneden, Cowan, & Gallino (2008) (83). All of the 
element abundances are normalized to the individual stellar log (e^;^) values (Eu is assumed 
to be an indicator of r-process contribution). For the n-capture elements with Z = 64 — 72, 
the RGB stellar abundance values strongly correlate with the solar r-process distribution. 
Similarly for the RHB stars, these abundances match well to the solar r-process pattern for 
most of the elements in the Z = 64 — 72 range. 



Figure [8] also displays the scaled, solar s-process abundance distribution (green, dotted 
line). The s-process predictions are also taken from Sneden, Cowan, & Gallino and the 
values are normalized to the solar log (esa) (Ba is considered to be an indicator of s-process 
contribution). As shown for the Z = 64 — 72 elements, there is virtually no agreement 
between the solar s-process pattern and either the RGB or the RHB stellar abundances. The 
s-process predictions compare well to the RGB abundances for only two elements: Ce and 
La. Thus, it follows that the nucleosynthesis of the heavy neutron capture elements in M15 
was dominated by the r-process. In addition, the abundance pattern for the light n-capture 
elements (Sr, Y, Zr) does not adhere to either a solar r-process or s-process distribution. 



5.1. Evidence for Additional Nucleosynthetic Mechanisms Beyond the 

Classical r- and s-process 

To further examine the anomalous light n-capture abundances in the M15 cluster. Fig- 
ures [9] and 10 are generated. For the stars of current effort and those from Otsuki et al. 



(2006) (jTTj), these two plots display the abundances of the n-capture elements (Sr, Y, Zr and 
La) as a function of the [Ba/H] and [Eu/H] ratios, respectively. Moreover, the abundance 
results from five select field stars which represent extremes in r-process or s-process enhance- 
ment are plotted (CS 22892-052 [Sneden et al. 2003, 2009(l8l ET])]; CS 22964-161 [Thompson 
et al. 2008](ll]); HD 115444 [Westin et al. 2000(j97D, Sneden et al. 2009(1811)]; HD 122563 
[Cowan et al. 2005(|23|, Lai et al. 2007([60l)]; HD 221170 [Ivans et al. 2006(i49j), Sneden et 



^B028, B412, and B584 are not included in the figure as these stars lack abundances for most of the 
elements in the specified Z range. 
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al. 2009(87)]). 

In Figure [9} an anti-correlative trend is seen for Sr, Y and Zr with Ba while no explicit 
correlative behavior is apparent for La. The correlation coefficient, r, is indicated in each 



panel. Likewise, La and Eu appear un-correlated in Figure 10 The [(Sr, Y, Zr)/Eu] ratios 
all exhibit anti-correlation with [Eu/H] in this figure. As shown, the elements Sr, Y, and Zr 
clearly demonstrate an anti-correlative relationship with both the markers of the s-process 
(Ba) and the r-process (Eu). 



Figures [9] and [TO] collectively imply that the production of the light neutron capture 
elements most likely did not transpire via the classical forms of the s-process or the r-process. 
This finding is not novel. The abundance survey of halo field stars by Travaglio et al. 
(2004) (l92il previously established the decoupled behavior of the light n-capture species to 
both Ba and Eu. Further, they postulated that an additional nucleosynthetic process was 
necessary for the production of these elements (Sr, Y , Zr) in metal-deficient regimes (coined 
the Lighter Element Primary Process; LEPP). 



The overabundances of Sr and Zr (see Figure |8j) could have been the result of a small 
s-process contribution to the M15 proto-cluster environment. To investigate this possibility, 
an abundance determination is performed for Pb, a definitive main s-process product. The 
upper panel of Figure [TT] illustrates the synthetic spectrum fits to the neutral Pb transition 
at a wavelength of 4057.8 A in the M15 giant, K462. An upper limit of log e{Pb) < -0.35 
can only be established for this star. For the remaining two RGB targets, upper limits were 
also determined and accordingly for all three, the average values of log e{Pb) < -0.4 and 
<[Pb/Eu]> < -0.15 were found. 



The lower panel of Figure 11 plots [Pb/Eu] as a function of [Eu/Fe] for the three M15 
RGB stars and the five, previously-employed halo field stars. In a recent paper, Roederer 
et al. (2010fl77jl) suggest that detections of Pb and enhanced [Pb/Eu] ratios should be strong 
indicators of main s-process nucleosynthesis. In turn, they contend that non-detections of 
Pb and depleted [Pb/Eu] ratios should signify the absence of nucleosynthetic input from 
the main component of the s-process (see their paper for further discussion). With the 
abundances of 161 low-metallicity stars ([Fe/H] < -1), Roederer et al. empirically determined 
a threshold value of [Pb/Eu] = -1-0.3 for minimum AGB contribution. As shown in the figure, 
the M15 giants lie below this threshold and accordingly, are likely devoid of main s-process 
input. Thus in the case of the M15 globular cluster, the light neutron capture elements 
presumably originated from an alternate nucleosynthetic process {e.g. u — p process, Frohlich 
et al. 2006 (j32!); high entropy winds, Farouqi et al. 2009(129])). 
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5.2. M15 Abundances in Relation to the Halo Field 



The upper panel of Figure 12 displays the evolution of the [Mg/Fe] abundance ratio 
with [Fe/H] for all M15 targets as well as for a sample of hundreds of field stars. For this 
figure, halo and disk star data have been taken from these surveys: Fulbright (2000) (|34j), 
Reddy et al. (2003) (jTED, Cayrel et al. (2004) (EQl), Cohen et al. (2004) (EI!) , Simmerer et 
al. (2004) (EnD, Barklem et al. (2005)®, Reddy, Lambert & Allende Prieto (2006) (75), 
Frangois et al. (2007) ([3lD, and Lai et al. (2008) ([59D. As shown, the scatter in the [Mg/Fe] 
abundance ratio is fairly small: A{[Mg/Fe\) m ax ~ 0.6 dex for all stars under conisderation 
and A([Mg/Fe])MAx ~ 0.1 dex for the M15 data set. In the metallicity regime below 
[Fe/H] < -1.1, the roughly consistent trend of [Mg/Fe] abundance ratio is due in part to 
the production history for these elements: magnesium originates from hydrostatic burning 
in massive stars while iron is manufactured by massive star, core-collapse SNe. If the short 
evolutionary lifetimes of these massive stars are taken into context with the abundance 
data, it would seem to indicate that the core-collapse SNe are rather ubiquitous events in 
the Galactic halo. Accordingly, the products that result from both stellar and explosive 
nucleosynthesis of massive stars should be well-mixed in the interstellar and intercluster 
medium. The apparent downward trend in the [Mg/Fe] ratio, in the metallicity region with 
[Fe/H]> -1.1, is due to nucleosynthetic input from Type la SNe, which produce much more 
iron in comparison to Type II events. 



In a similar vein, the lower panel of Figure 12 plots [Eu/Fe] as a function of [Fe/H 



and demonstrates that as the metallicity decreases, the spread in the [Eu/Fe] abundance 



ratio increases enormousl}- By contrast, the scatter in the M15 [Eu/Fe] ratios is large and 
comparable to the spread of the halo field at that metallicity. Specifically in the metallicity 
interval -2.7 < [Fe/H] < -2.2, the scatter in the [Eu/Fe] ratio is found to be cr = ±0.27 
for the nine stars of the M15 sample and similarly for the 23 halo giants, the associated 
scatter is cr = ±0.33. This variation in the relative europium abundance ratio (as first 
detected by Gilroy et al. 1988( 13^ and later confirmed by others, e.g. Burris et al. 2000(14), 
Barklem et al. (2005) (IE])) indicates an inhomogeneous production history for Eu and other 
corresponding r-process elements. These elements likely originate from lower mass SNe and 
their production is not correlated with that of the alpha elements (Cowan & Thielemann 
2004(24p). Furthermore, it seems that nucleosynthetic events which generated the r-process 
elements were rare occurrences in the early Galaxy. As a consequence, these elements were 
not well-mixed in the interstellar and intercluster medium (Sneden et al. 2009(87)). Note 



■'Though the data sample of Figure 



12 



is compilation of several sources, the scatter in the [Mg/Fe] and 



[Eu/Fe] ratios duplicates that found by such large scale surveys as, e.g. Barklem et al. (2005) 
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that r-process enhancement seems to be a common feature of all globular clusters (e.g. 
Gratton, Sneden, & Carretta 2004(36)). On the other hand, the scatter in select r-process 
element abundances, as found in M15, is not. 

6. SUMMARY 

A novel effort was undertaken to perform a homogenous abundance determination in 
both the RGB and RHB members of the M15 globular cluster. The current investigation 
employed improved atomic data, stellar model atmospheres, and radiative transfer code. A 
resolute offset in the iron abundance between the RGB and RHB stars on the order of 0.1 
dex was measured. Notwithstanding, the major findings of the analysis for both the RGB 
and RHB stellar groups include: a definitive r-process enhancement; a significant spread in 
the abundances of the neutron capture species (which appears to be astrophysical in nature); 
and, an anti-correlation of light n-capture element abundance behavior with both barium 
([Ba/H]) and europium ([Eu/H]). Accordingly, the last set of findings may offer proof of 
the operation of a LEPP-type mechanism within M15. To determine if these abundance 
behaviors are generally indicative of very metal-poor globular clusters, a comprehensive 
examination of the chemical composition of the analogous M92 cluster should be undertaken 
([Fe/H] ~ -2.3, Harris et al. 1996(H?I); the literature contains relatively little information 
with regard to the n-capture abundances for this cluster). 

To date, the presence of multiple stellar generations within the globular cluster M15 has 
not been irrefutably established. In a series of papers, Carretta et al. (2009a, 2009c, 2010(1T^ 
[TB| [T7|) ) offered compelling proof in the detection of light element anti-correlative behavior 
(Na-0) in numerous members of the M15 RGB. Lardo et al. (2011 (l6Tll ) did find a statistically 
significant spread in the SDSS photometric color index of u — g, but yet they were not able 
to demonstrate a clear and unambiguous correlation of {u — g) with the Na abundances in 
the RGB of the M15 cluster (which would have provided further evidence). To wit, recent 
investigations of M15 have revealed several atypical features including: probable detection 
of an intermediate mass black hole (van der Marel et al. 2002(93); though the result is under 
some dispute); observation of an intracluster medium (Evans et al. 2003 (l28|l ): detection of 
mass loss (Meszaros et al. 2008, 2009(1^7} IS^ ): identification of extreme horizontal branch 
and blue hook stars (Haurberg et al. 2010 (H5]1 ): and, observation of an extended tidal tail 
(Chun et al. 2010fi21i)). It would be worthwhile to examine these peculiar aspects of the 
globular cluster in relation to the abundance results of M15. Further scrutiny is warranted in 
order to understand the star formation history and mixing timescale of the M15 protocluster 
environment. 
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A. Synopsis 

The essential approacli to tlie solution of radiative transfer in the MOOG program 
has been altered with the employment of short characteristics and the application of an 
accelerated lambda iteration (ALI) scheme. Original development of the short characteris- 
tics (SC) methodology in the context of radiative transfer was done by Mihalas, Auer, & 
Mihalas (1978( 1681) ). Improvement of the SC approach in the explicit specification of the 
source function (at all grid points) was made by Olson & Kunasz (1987) (170ll and Kunasz & 
Auer (1988 (13Hll ). As a supplemental reference, the current version of MOOG draws upon 
the concise treatment of Auer & Paletou (1994([7])). The implementation of the ALI tech- 
nique within the framework of radiative transfer and stellar atmospheres was first done by 
Werner (1986f l96|l ) and subsequently refined by both e.g.^ Rybicki & Hummer (1991(78!)) 
and Hubeny (1992 (Il8|) ). The main SC and ALI prescription followed by the MOOG code 
is that from Koesterke et al. (2002 (I52p and references therein). Since the text below is a 
general description and pertains to the specific coding in MOOG, the reader should consult 
the aforementioned references as they contain significantly more information. 

The primary modification to MOOG is the creation and incorporation of four new 
subroutines. The names and purposes of the new subroutines are as follows: AngWeight . f , 
which determines the Gaussian weights and integration points; Sourcefunc_scat_cont . f , 
which incorporates both a scattering and an absorption component to compute the source 
function and resultant fiux for the continuum; Sourcef unc_scat_line . f , which incorporates 
both a scattering and an absorption component to compute the source function and resultant 
flux for the line; and, Cdcalc_JS.f, which calculates the final line depth via the emergent 
continuum and line fluxes. To accommodate these additions, several key subroutines were 
also revised. Further details and the publicly-available MOOG code may be found at the 
website: http://verdi.as.utexas.edu/moog.html, Note that MOOG still retains the 
capacity to operate in -pure absorption mode with the source function set simply to S = eB. 



B. Contributions to the Continuous Opacity 

To remind the reader, in the visible spectral range, the two principal sources of opacity 
in stellar atmospheres are the bound-free absorption from the negative hydrogen ion (H~) 
and Rayleigh scattering from neutral atomic hydrogen. The standard expression {e.g. Gray 
1976 (13S|) ) for the B.~ bf absorption coefficient is 

n ~ 4.1458 X 10-^''aBFPeQ''^^W''-^^^'^ (Bl) 
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where a^F is the bound-free atomic absorption coefficient (which has frequency dependence), 
Pe is the electron pressure, and G = 5040/T (note that Eg. |Bl| is per neutral hydrogen atom). 

With regard to Rayleigh scattering, the scattering cross-section of radiation with angular 
frequency u incident upon a neutral H atom is given by the Kramers-Heisenberg formula in 
terms of atomic units as 



a(w) ^ 



(B2) 



where ctt is the Thompson scattering cross-section and Ui is the angular frequency corre- 
sponding to the Lyman limit. Numerical calculations of the Ai coefficients and the generation 
of an exact expression for Eq. |B2| have been done by Dalgarno & Williams (1962 (l26jl ) and 
more recently by Lee & Kim (2004(63)). 

The H^BF and Rayleigh scattering opacity contributions depend on temperature and 
metallicity (and to some extent, on the surface gravity). Rayleigh scattering also has a 
dependence, and as a consequence, it greatly influences blue wavelength transitions. For the 
majority of stars (such as dwarfs and sub-giants), B.~ bf is the dominant opacity source in the 
visible spectral regime. However for low temperature, low metallicity giants, the Rayleigh 
scattering contribution becomes comparable to and even exceeds that from B.^ bf in the 
ultraviolet and blue visible wavelength regions. Therefore, to accurately determine the line 
intensity with the correct amount of flux and opacity contribution for all stellar types and 
over a wide spectral range, isotropic, coherent scattering must be considered. 



C. Form of the Radiative Transfer Equation and Implementation of the ALI 

Scheme 

The source function is then written as: S = {1 — e)J + eB (where e is the thermal 
coupling parameter, J is the mean intensity, and B is the Planck function). To commence 
with the formal solution of the radiation transfer equation, a fundamental assumption is 
made in that the source function is specified completely in terms of optical depth. After 
some mathematical manipulation, the radiative transfer equation becomes 



,d^ 
dr^ 



J-S 



(CI) 



where /i is the directional cosine and r is the optical depth. Eq. CI is an integro-difFerential 
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equation (and subject to boundary conditions). To obtain the numerical solution of Eq. CI 



a discretization in angle and optical depth is necessary. As a consequence, the solution is 
simplified and a Gaussian quadrature summation is done instead of an integration. Though 
it is eventually possible to evaluate Eq. CI in a single step, the use of an iterative method 
to arrive at a solution is preferred as it is computationally faster than a straightforward 
approach. For the MOOG program, the ALI technique is employed with the application of a 
full acceleration. In the context of ALI scheme, the transfer equation takes the form of J = 
A[S], where A represents the matrix operator. Through the concept of preconditioning, ALI 
allows for the efficient, iterative solution of a (potentially) large system of linear equations. 
The main steps of the iterative cycle are the evaluation of J = A[S'], the computation of the 
AS quantity, and the corresponding adjustment to the source function. 



D. Short Characteristic Solution of Line Transfer 



From a general viewpoint, radiative transfer can be thought of as the propagation of 
photons along a ray on a two-dimensional grid. Note that the number of rays corresponds to 
the number of quadrature angles. Determination of radiation along a ray is done periodically 
at ray segments, or short characteristics. Essentially, SC start at a grid point and proceed 
along the ray until a cell boundary is met. At these cell boundaries, the intensity, opacity, 
and source function values are established. Then with the knowledge of the cell boundary 
intensities, the intensity at other non-grid points can be calculated. 

Specifically with regard to MOOG, the intensity determination is a function of depths 
(i) and angles (j). It is performed for both an inward (i — l) and an outward (i + 1) direction. 
Along the characteristic, the opacity quantity is assumed to be a linear function. In effect, 
the intensity for the ray can be expressed as 



/ = /e-^"« + / S{T)e-^dT (Dl) 



The optical depth step At can be thought of as the path integral of the opacity 
along the characteristic (the entire, involved definition of the At quantity is found in the 



Sourcef unc_scat_* subroutines). Now, the evaluation of Eq. Dl requires the interpolation 
of the source function. A linear interpolation is sufficient to satisfy the various boundary 
conditions. Interpolation over [S{i), S{i =F 1)] then entails 



S{T)e^^dT = S{i)wi{i) + S{i T l)w2(i) 



(D2) 



The weights are given by the relations 



^o(^) = (e-^^« - 1)/Ar(z) (D3) 
w^{{) = 1 + Wo{i) (D4) 
= -e-^-(') - wo{i) (D5) 

These weights are found by recursion. The use of hnear interpolation does not generate 
significant error (as normally would occur) due to the optically-thin nature of the boundary 
layer. The expression for the mean intensity, J, subsequently becomes 



where the Wcau are the Gaussian quadrature weights (these are distinct from the weights of 



Eq. D2). The summation over all depth points and angles(/rays) is necessarily performed. 



The SC formal solution of the transfer equation then proceeds in an iterative manner. 

Ongoing and future improvements to the MOOG code include: the incorporation of 
the Lee & Kim (2004 (IBHI) ) formulation for Rayleigh scattering (for atomic hydrogen), the 
employment of a further discretization with regard to frequency, and the implementation of 
spherically-symmetric geometry in the solution of radiative transfer. 
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Fig. 1. — A color-magnitude diagram for the globular cluster M15. For cluster members 
(indicated by small black dots), V magnitudes are taken Buonanno et al. (1983) while the 
K magnitudes are obtained from the 2MASS database (Skrutskie et al. 2006). Note that 
stars with both B and V magnitudes from Buonanno et al. are shown in the plot. The 
program stars are signified by large, red circles. Also displayed are the isochrone data from 
both Marigo et al. (2008; labeled PM2008; shown in green) and Dotter et al. (2008; labeled 
AD2008; shown in blue). 
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Fig. 2. — A small wavelength region shown for all program stars. For the purposes of display, 
the relative fluxes of target stars (other than K341) have been shifted by additive constants. 
Vertical dotted lines denote the spectral features of La II at a wavelength of 4123.22 A and 
of Eu II at 4129.72 A. Note that these transitions appear to be more pronounced in the 
spectra of the lower temperature RGB targets. 
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Fig. 3. — Binned abundances for the elements with Z = 6,8,12 — 14,20 — 21 displayed 
as quartile boxplots. For each species, the boxplot shows the median, the upper and lower 
quartiles, and the extremes for both the RGB (denoted in red) and RHB (signified in blue) 
samples. Outliers, data points with values greater than 1.25 times the median, are indicated 
by star symbols. Note that there are no RHB points for C I and O I. Note that a depletion 
in the relative carbon abundance with respect to solar is found. 
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Fig. 4. — Binned abundances for the elements with 22 < Z < 27 displayed as quartile 
boxplots. For each species, the boxplot shows the median, the upper and lower quartiles, 
and the extremes for both the RGB (denoted in red) and RHB (signified in blue) samples. 
Note that no RHB data exist for V I. Also, the Mn I results have been set aside for the RHB 
stars (consult the text for further information). The most consistent agreement between 
neutral and first-ionized species is found for Ti in all M15 stars. 
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Fig. 5. — Binned abundances for the elements with Z = 29,30,38 — 40,56 — 59 displayed 
as quartile boxplots. For each species, the boxplot shows the median, the upper and lower 
quartiles, and the extremes for both the RGB (denoted in red) and RHB (signified in blue) 
samples. Outliers, data points with values greater than 1.25 times the median, are indicated 
by star symbols. Note that there are no RHB data points for Cu I and also, only one RHB 
abundance value for Zn I. 



- 40 - 



■ii -1=1=1- jD 

1 I □ >H 



B 



o _ 



Q 



X3 



r-|— I 73 
I I CJ 



3 



T3 _ 

z 



c 

e 

3 



If5 O lO O lO O 

cv? in o w in 

^ ^ ^ cvi c\j cvi 

I I I I 



[H/lU8Uiaia] 



Fig. 6. — Binned abundances for the elements with 60 < Z < 70 displayed as quartile 
boxplots. For each species, the boxplot shows the median, the upper and lower quartiles, 
and the extremes for both the RGB (denoted in red) and RHB (signified in blue) samples. 
Outliers, data points with values greater than 1.25 times the median, are indicated by star 
symbols. Note that only one RHB abundance value was found for the elements Nd II, Sm 
II, and Tm II. In the case of Tb II, no RHB measurements were possible. 
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Fig. 7. — Plots of [(Ba, La)/H] as a function of [Eu/H] for the M15 stars of the current 
study as well as those from the original Sneden et al. (1997) paper. Also shown are the 
associated error bars for the abundances (these are the standard deviation values given in 
Table |4]). Note that the Eu and La abundance determinations for the 1997 stellar sample 
are based upon one transition only (and consequently, no error bars are plotted). Contrary 
to Sneden et al., no clear evidence of a binary distribution in the abundances is detected. 
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Fig. 8. — Comparison of the n-capture abundances of six M15 targets to the Solar System r- 
and s-only abundance distributions. The abundances of the three RGB stars (K341, K462, 
and K583) are signified by red symbols while the abundances of the three RHB stars (B009, 
B224, and B262) are designated by blue symbols. The solid black line denotes the r-process 
only abundance pattern, which is scaled to the solar log e{Eu), and the dotted green line 
indicates the s-process only abundance pattern, which is scaled to the solar log e{Ba) (all 
predictions are taken from Sneden, Cowan, & Gallino 2008). The heavy element abundances 
of the M15 stars compare well to the r-process predictions, but not those of the s-process. 
Neither abundance distribution consistently matches the stellar abundances for the elements 
Sr, Y, and Zr. 
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Fig. 9. — Plot of the relative abundances of the light n-capture elements Sr, Y, and Zr as 
function of the traditional s-process indicator, Ba. All M15 stars from the current effort 
are shown. Also displayed are the abundance results from Otsuki et al. (2006) as well as 
those for a few halo field stars (taken from various literature references; see text for further 
details). The correlation coefficient, r, is exhbited in each of the panels and as seen, Sr, Y , 
and Zr demonstrate a clear anti-correlative behavior with Ba while La does not. 
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Fig. 10. — Plot of the relative abundances of the light n-capture elements Sr, Y, and Zr as 
function of the traditional r-process indicator, Eu. All M15 stars from the current effort 
are shown. Also displayed are the abundance results from Otsuki et al. (2006) as well as 
those for a few halo field stars (taken from various literature references; see text for further 
details). The correlation coefficient, r, is exhbited in each of the panels and as seen, Sr, Y , 
and Zr demonstrate a clear anti-correlative behavior with Eu while La does not. 
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Fig. 11. — Upper panel: Comparison of synthetic and observed spectra for the Pb I transition 
at 4057.8 A in the M15 K462 giant. Four incremental changes in abundance are shown for 
the specified Pb feature. It is only possible to establish an upper limit of approximately log e 
(Pb)~ -0.35 for this star. A CH contaminant is present in the blue wing of the Pb transition 
and accounts for a definitive portion of the signal. Accordingly, it appears that the s-process 
element Pb is nominally detected in K462. Lower Panel: Plot of [Pb/Eu] as a function of 
[Eu/Fe] for three M15 giants and 5 halo field stars. The empirically determined threshold 
ratio, which indicates the occurrence of AGB enrichment, is shown by the short-dashed line 
at [Pb/Eu] > -1-0.3 (this value is taken from Roederer et al. 2010; see the text for further 
details). Notice that all of the M15 stars fall below this line (and correspondingly, should 
lack main s-process material). 
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Fig. 12. — Comparison of the distribution [(Mg, Eu)/Fe] relative abundances as a function 
of metallicity for M15 targets (denoted by the blue and red circles) as well as for halo and 
disk stars (signified by the green, filled circles). Field star data have been taken from large 
sample surveys (see the text for further information). As per convention, the dotted lines 
represent the Solar abundance ratios. In in the lower panel, the two dashed black lines 
indicate the approximate ranges of the [Eu/Fe] data. The spread in the [Mg/Fe] values is 
smaller for M15 than that for the field. Yet for [Eu/Fe], the M15 scatter is comparable to 
that of the field (for that particular metallicity). 
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Table 4. Average Relative Abundance Data for M15 Program Stars 





K341 


K462 


K583 


B009 


B028 


B224 


B262 


B412 


B584 


<[C7/Fe7]> 


-0.34 


-0.08 


-0.40 


-0.10 












a (loge(C,)) 




















No. of Lines 


1 


1 
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1 

















<[0//Fe7]> 


0.85 


0.78 


0.62 














o (loge(0/)) 


0.01 


0.09 
















No. of Lines 


2 


2 


1 




















<[Na//Fe/]> 


-0.18 


-0.11 


1.11 


1.60 


0.46 


1.18 


1.21 


0.25 


0.28 


u (loge(Naj)) 




0.13 


0.00 


0.01 


0.12 


0.09 


0.04 


0.01 


0.08 


No. of Lines 
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2 
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<[Mg7/Fcj]> 


0.60 


0.45 


0.48 


-0.01 


0.56 


0.56 


0.44 


0.57 


0.54 


o (loge(Mg/)) 


0.17 


0.13 


0.07 


0.09 


0.15 


0.07 


0.20 


0.23 


0.16 


No. of Lines 


3 
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4 
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2 


1 


2 


2 
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<[Alj/Fc/]> 


0.23 


0.29 


0.64 


0.34 


-0.60 


-0.27 


-0.26 


-0.52 


-0.52 


o (loge(Al,)) 


0.18 


0.17 


0.16 


0.11 


0.14 


0.08 


0.22 
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No. of Lines 
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<[Si//Fe/]> 
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a (loge(Ca/)) 
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No. of Lines 
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<[Sc//Fe/]> 






-0.76 
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No. of Lines 








1 




















<[Sc///Fe,,]> 


0.16 


0.04 


0.12 


0.37 


0.43 


0.37 


0.16 


0.39 


0.32 


a (loge(Sc,,)) 


0.10 


0.12 


0.10 


0.16 


0.06 


0.08 


0.09 


0.20 




No. of Lines 


2 


2 


2 


2 


3 


3 


4 


2 


1 


<[Ti//Fe/]> 


0.11 


0.19 


0.07 


0.58 


0.60 


0.54 


0.32 


0.85 


0.77 


o (loge(Ti/)) 


0.12 


0.13 


0.11 


0.16 


0.16 


0.15 


0.15 


0.08 


0.14 


No. of Lines 


15 


12 


13 


2 


4 


4 


6 


2 


2 


<[Tijj/Fo,,]> 


0.47 


0.44 


0.47 


0.38 


0.48 


0.40 


0.31 


0.52 


0.50 


CT (loge(Ti//)) 


0.11 


0.13 


0.14 


0.10 


0.10 


0.10 


0.07 


0.07 


0.06 


No. of Lines 


14 


8 


16 


9 


7 


9 


11 


8 


7 


<[V7/Fe7]> 


-0.26 


-0.11 


-0.18 














o (loge(V/)) 


0.07 


0.19 


0.12 














No. of Lines 


2 


3 


3 
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Table 4 — Continued 





K341 


K462 


K583 


B009 


B028 


B224 


B262 


B412 


B584 


<[V///Fe//]> 


0.27 


0.43 


0.11 


0.23 


0.41 


0.31 


0.04 




0.44 


a (loge(V//)) 


0.01 


0.21 


0.21 




0.05 




0.04 




0.03 


No. of Lines 


2 


2 


2 


1 


2 


1 


2 





2 


<[Cr7/Fej]> 


-0.36 


-0.22 


-0.38 


-0.21 


-0.24 


-0.21 


-0.22 


-0.07 


-0.19 


a (loge(Crj)) 


0.05 


0.12 


0.07 


0.15 


0.05 


0.12 


0.21 


0.10 


0.09 


No. of Lines 


5 


6 


5 


4 


3 


3 


8 


3 


2 


<[Cr,,/Fe,/]> 


-0.06 


-0.02 


-0.03 


0.33 


0.37 


0.24 


0.10 


0.43 


0.38 


a (loge(Crj,)) 


0.14 


0.18 


0.06 


0.01 


0.21 




0.10 


0.06 


0.07 


No. of Lines 


4 


3 


3 


2 


3 


1 


3 


2 


2 


<[Mn//Fe/]> 


-0.34 


-0.45 


-0.42 


-0.70 


-0.74 


-0.80 


-0.78 






a (loge(Mn7)) 


0.02 


0.04 
















No. of Lines 


2 


2 


1 


1 


1 


1 


1 








<[Fe//H]> 


-2.53 


-2.54 


-2.57 


-2.69 


-2.65 


-2.74 


-2.77 


-2.62 


-2.63 


a (loge(Fe7)) 


0.14 


0.12 


0.14 


0.12 


0.12 


0.11 


0.14 


0.12 


0.15 


No. of Lines 


59 


57 


63 


98 


72 


72 


100 


45 


60 


<[Fcjj/H]> 


-2.51 


-2.51 


-2.57 


-2.65 


-2.60 


-2.66 


-2.73 


-2.59 


-2.59 


a (loge(Fc„)) 


0.13 


0.11 


0.10 


0.09 


0.07 


0.13 


0.13 


0.07 


0.13 


No. of Lines 


8 


9 


10 


12 


8 


9 


11 


6 


9 


<[Co//Fe/]> 


0.07 


0.25 


0.35 


0.45 


0.18 


0.18 


0.04 




0.22 


a (loge(Co/)) 


0.21 


0.18 


0.20 


0.32 












No. of Lines 


2 


2 


2 


2 


1 


1 


1 





1 


<[Ni//Fe/]> 








-0.23 












a (loge(Ni,)) 








0.11 












No. of Lines 











2 

















<[Cu//Fe/]> 


-0.92 


-0.88 


-0.92 














a (loge(Cu,)) 




















No. of Lines 


1 


1 


1 




















<[Zn//Fc/]> 


0.02 


-0.01 


0.01 






0.49 


0.27 






a (loge(Zn/)) 


0.07 


0.07 










0.06 






No. of Lines 


2 


2 


2 








1 


2 








<[Sr„/Fcjj]> 


-0.01 


0.28 


-0.18 


0.16 


0.35 


0.43 


0.06 


0.38 


0.46 


cr (loge(Si7/)) 


0.00 


0.04 


0.11 


0.04 


0.11 


0.14 


0.11 


0.14 


0.11 


No. of Lines 


2 


2 


2 


2 


2 


2 


2 


2 


2 


<[Y77/Fejj]> 


-0.13 


-0.04 


-0.23 


0.19 


0.02 


0.16 


-0.04 


0.18 


0.08 


a (loge(Y,,)) 


0.04 


0.04 


0.04 


0.06 


0.11 


0.09 


0.09 


0.04 


0.03 


No. of Lines 


7 


7 


6 


6 


2 


4 


10 


2 


3 
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Table 4 — Continued 





K341 


K462 


K583 


B009 


B028 


B224 


B262 


B412 


B584 


<[Zr///Fc//]> 


0.21 


0.27 


0.14 


0.53 


0.35 


0.38 


0.35 




0.71 


a (loge(Zr„)) 


0.06 


0.12 


0.17 


0.16 




0.14 


0.02 




0.00 


No. of Lines 


4 


4 


4 


4 


1 


2 


5 





2 


<[Ba77/Fe7j]> 


0.20 


0.23 


0.01 


0.40 


-0.14 


0.21 


0.43 


0.14 


-0.03 


a (loge(Ba/j)) 


0.03 


0.10 


0.07 


0.08 


0.03 


0.05 


0.10 


0.07 


0.06 


No. of Lines 


4 


4 


4 


5 


3 


5 


6 


3 


3 


<[La///Fe//]> 


0.16 


0.42 


-0.04 


0.60 


-0.01 


0.58 


0.38 




0.60 


a (loge(La77)) 


0.08 


0.09 


0.09 


0.07 




0.10 


0.08 




0.14 


No. of Lines 


9 


9 


9 


8 




6 


11 


1 


2 


<[Ce77/Fe77]> 


0.20 


0.33 


-0.17 


0.51 






0.37 






a (loge(Ce77)) 


0.15 


0.14 


0.13 


0.13 






0.14 






No. of Lines 


25 


22 


11 


10 








13 








<[Prj7/FC77]> 


0.42 


0.63 


0.28 


0.73 




0.81 


0.69 






a (loge(Pr77)) 


0.04 


0.03 


0.12 


0.00 






0.11 






No. of Lines 


5 


5 


5 


2 





1 


2 








<[Nd77/Fe77]> 


0.47 


0.68 


0.23 


0.58 












a (loge(Nd77)) 


0.10 


0.11 


0.14 


0.10 












No. of Lines 


27 


25 


22 


2 

















<[Sm77/Fe77]> 


0.50 


0.77 


0.29 


0.61 












a (loge(Sm77)) 


0.08 


0.08 


0.10 


0.13 












No. of Lines 


20 


20 


15 


3 

















<[EU77/Fe77]> 


0.50 


0.83 


0.26 


0.90 


0.72 


0.85 


0.69 


1.18 


0.96 


a (loge(Eu77)) 


0.12 


0.11 


0.11 


0.03 


0.04 


0.09 


0.07 


0.07 


0.18 


No. of Lines 


7 


6 


7 


5 


2 


5 


5 


2 


2 


<[Gd77/Fe77]> 


0.55 


0.77 


0.25 


0.69 






0.67 






a (loge(Gd77)) 


0.08 


0.07 


0.03 














No. of Lines 


8 


9 


3 


1 








1 








<[Tb///FC77]> 


0.56 


0.67 


0.42 








0.96 






a (loge(Tb77)) 




















No. of Lines 


1 


1 


1 











1 








<[Dy77/FC77]> 


0.38 


0.58 


0.36 


1.13 






0.75 






cr (loge(Dy//)) 


0.15 


0.03 


0.07 


0.08 






0.11 






No. of Lines 


4 


2 


4 


3 








3 








<[H077/Fe77]> 


0.45 


0.67 


0.26 


1.09 




1.35 


0.81 






a (loge(Ho77)) 


0.05 


0.06 


0.10 


0.10 






0.05 






No. of Lines 


3 


4 


3 


3 





1 


4 
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Table 4 — Continued 





K341 


K462 


K583 


B009 


B028 


B224 


B262 


B412 


B584 


<[Er///Fe//]> 


0.48 


0.67 


0.18 


0.88 




0.73 


0.58 






a (loge(Er//)) 


0.15 


0.15 


0.06 


0.03 




0.14 


0.05 






No. of Lines 


4 


4 


3 


4 





2 


3 








<[Tm///Fe//]> 


0.40 


0.66 


0.26 








0.85 






£7 (loge(Tni7/)) 


0.09 


0.20 


0.07 














No. of Lines 


3 


3 


2 











1 








<[Yb/j/Fc//]> 


0.23 


0.64 


0.09 


0.93 


1.29 


0.71 




0.89 


0.71 


a (loge(Yb77)) 




















No. of Lines 


1 


1 


1 


1 


1 


1 





1 


1 


<[Hf,,/Fc„]> 


0.21 


0.52 


-0.08 














a (loge(Hf//)) 




















No. of Lines 


1 


1 






















<[Os7/Fe7]> 














2.00 






a (loge(Osz)) 














0.35 






No. of Lines 




















2 








<[Ir//Fe/]> 














1.35 






£7 (loge(Ir7)) 




















No. of Lines 




















1 








<[Pb//Fe/]> 


<0.44 


< 0.35 


< 0.38 














a (loge(Pb7)) 




















No. of Lines 


1 


1 


1 




















<[Th/7/Fe77]> 


0.46 


0.67 


0.32 














a (loge(Th7/)) 




















No. of Lines 


1 


1 


1 
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Table 5. Solar Photospheric Abundances of the Elements Analyzed in the Current Study 



Element Z log{e ei) Reference Methodology 



"LJ 

ri 


1 
i 


12.00 


M / A 








D 


8.43 it 


0.05 


Aspiuna et ai. zuuy 


Q "Pi 

o-U 






Q 
O 


8.71 i 


0.05 


ocott et al. zuuy 


o-U 




IN a 


1 1 


6.1 / ± 


0.04 


/\spiuna ci ai. zuuo 


o-U 




Mg 


1 o 


7.53 ± 


0.09 


Aspiuna et al. zuuo 


Q T\ 
O-U 


T 'FT? 


A 1 




Q*? _1_ 


0.U6 


Aspiuna et al. zuuo 


o-L) 


Lil III 


ol 


1 A 

141: 


7.51 ± 


0.U4 


Aspiuna et ai. ^uuo 


o-U 


T TTT 

Lil III 




on 


6.31 ± 


0.04 


Aspiuna et al. zuuo 


O-U 


T "FT? 
Lil III 


GC 




3.15 ± 


0.04 


Aspiuna ci ai. zuuy, ^.jrcvcssc ci ai. zuiu 


O-U 


T "FF 


1 1 


oo 


4.95 ± 


0.05 


Am-.ln»^..l ol onno. r^titi^ ot- ol omn 
Aspiuna ct al. zUUy, Lrrevesse et al. zUiU 


O-U 




V 


zo 


3.93 ± 


0.08 


AnT^ln»^..l ol onno. f~^\-f^^m^^n ol 0010 

Aspiuna ct al. zUUy, Lrrevesse et al. zUiU 


O-U 




or 


Z4 


5.64 ± 


0.02 


ooDecK et al. zuu/ 


1-U 




Ayr 1-1 


zo 


5.43 ± 


0.04 


AriT^lim^ ni- ol 0000. me^e^n o+ ol 0010 

Aspiuna et al. zuuy; tjievesse et ai. zuiU 


O-U 


T "FT? 
Lil III 


re 


ZD 


7.52 ± 


0.08 


C!__ _4. ol 1 001 

oneaen et ai. lyyi 


1 F\ 
i-U 


T "FT? 




07 
Z f 


4.99 ± 


0.08 


A o i^il 1 I Ti ..1 o4- ol OOOO. /^T.OTroooo o4- ol 0010 

Aspiuna et ai. zuuy; orevesse et ai. zuiU 


O-U 


T "FT? 


IV i 
iNl 


zo 


6.17 ± 


0.02 


ocott ct ai, zuuy 


O-U 


T "FTT 
J_j1 Hj 




zy 


4.19 ± 


0.04 


A&piuna ct ai. zuuy, i^icvcssc ct ai. zuiu 


Q Fl 
O-U 


J_j i r!j 


ZjII 


oU 


4.60 ± 


0.03 


oieinoiit v.joaciroia lyou 


1 Fl 
1-U 


T 'FF 


Sr 


38 


2.92 ± 


0.05 


Barklem & O'Mara 2000 


l-D 


LTE 


Y 


39 


2.21 ± 


0.05 


Asplund et al. 2009; Grevesse et al. 2010 


3-D 


LTE 


Zr 


40 


2.58 ± 


0.04 


Asplund ct al. 2009; Grevesse ct al. 2010 


3-D 


LTE 


Ba 


56 


2.17 ± 


0.09 


Asplund ct al. 2009; Grevesse ct al. 2010 


3-D 


LTE 


La 


57 


1.13 ± 


0.03 


Lawlcr ct al. 2001 


l-D 


LTE 


Ce 


58 


1.61 ± 


0.01 


Lawler et al.2009 


l-D 


LTE 


Pr 


59 


0.76 ± 


0.02 


Sneden et al. 2009 


l-D 


LTE 


Nd 


60 


1.45 ± 


0.01 


Den Hartog et al. 2003 


l-D 


LTE 


Sm 


62 


1.00 ± 


0.01 


Lawler et al. 2006 


l-D 


LTE 


Eu 


63 


0.52 ± 


0.01 


Lawlcr ct al. 2001 


l-D 


LTE 


Gd 


64 


1.11 ± 


0.01 


Den Hartog ct al. 2006 


l-D 


LTE 


Tb 


65 


0.28 ± 0.3 


Lawlcr ct al. 2001 


l-D 


LTE 


Dy 


66 


1.13 ± 


0.02 


Sneden et al. 2009 


l-D 


LTE 


Ho 


67 


0.51 ± 


0.10 


Lawler et al. 2004 


l-D 


LTE 


Er 


68 


0.96 ± 


0.03 


Lawler et al. 2008 


l-D 


LTE 


Tm 


69 


0.14 ± 


0.02 


Sneden et al. 2009 


l-D 


LTE 


Yb 


70 


0.86 ± 


0.10 


Sneden ct al. 2009 


l-D 


LTE 


Hf 


72 


0.88 ± 


0.08 


Lawlcr ct al. 2007 


l-D 


LTE 


Os 


76 


1.25 ± 


0.11 


Quinet et al. 2006 


l-D 


LTE 


Ir 


77 


1.38 ± 


0.05 


Youssef & Khalil 1988 


l-D 


LTE 


Pb 


82 


1.75 ± 


0.10 


Asplund et al. 2009; Grevesse et al. 2010 


3-D 


LTE 


Th 


90 


0.02 ± 


0.10 


Asplund et al. 2009; Grevesse et al. 2010 


3-D 


LTE 



